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We have calculated surface energies and surface magnetic order of various low-indexed surfaces of 
monoatomic Fe, Co, and Pt, and binary, ordered FePt, CoPt, and MnPt using density functional theory. Our 
results for the binary systems indicate that elemental, Pt-covered surfaces are preferred over Fe- and Co-covered 
and mixed surfaces of the same orientation. The lowest energy orientation for mixed surfaces is the highly 
coordinated (111) surface. We find Pt-covered (111) surfaces, which can be realized in the Lli structure only, 
to be lower in energy by about 400 meV/atom compared to the mixed LIq (111) surface. We conclude that in 
small nanoparticles this low surface energy can stabilize the Lli structure, which is suppressed in bulk alloys. 
From the interplay of surface and bulk energies, equilibrium shapes of single-crystalline ordered nanoparticles 
and crossover sizes between the different orderings can be estimated. 



I. INTRODUCTION 



During the last two decades, an exponential increase of 
the magnetic data storage areal density has been achieved. 
Thus, in order to continue with this trend, a constant fur- 
ther miniaturization of the bit size is required. Promising 
candidates for future ultra-high density storage media are 
Llo-ordered FePt or CoPt nanoparticles due to their extraor- 
dinary high magnetocristalline anisotropy in the bulk phase 
(FePt: Ku = 7•10^erg/cm^ CoPt: = 4.9-10^erg/cm3)^ 
The Llo-lattice structure is characterized by a tetragonal 
distortion of a few percent along the c-axis accompanied by 
an alternating stacking of elemental layers along the [001] 
direction (cf. Fig. [Hi. The intriguing properties of LIq FePt 
and CoPt alloys and nano-composites have been subject to 
numerous experimental and theoretical studies, e.g., see Refs. 
I llellTlliMlolfTlin^^ references therein. 

A subtle interplay between surface energies and internal inter- 
face energies determines the equilibrium shape of nanoparti- 
cles. One major obstacle in producing LIq FePt nanoparticles 
is the occurrence of multiple twinning i i^i i^i i "^i i 1 9.20,2 1 
Multiply-twinned nanoparticles such as icosahedra or decahe- 
dra do not exhibit high uniaxial magnetocrystalline anisotropy 
energy due to the different crystallographic orientation of the 
individually ordered twins. Multiple twinning appears if the 
energy gain due to low surface energies exceeds the energy 
needed for the creation of twin boundaries. The hierarchy of 
surface energies is thus one important function determining 
the equilibrium shape of small nanoparticles^i Since surface 
energies are particularly difficult to measure in experiment, 
their theoretical calculation is an important task. 
Apart from the LIq phase, also in the less common Lli 
structure a high uniaxial magnetic anisotropy is reported i^^'^** 
In the Lli structure, alternating fee Cu and Pt layers are 
stacked along the [111] direction, similar to the LIq structure, 
which consists of alternating (001) planes (cf. Fig. [T]). In 
contrast to the LIq structure, the Lli phase is only stable for 
bulk materials in the metallic CuPt alloy . Very recently, 
Lli type CoPt ordered films with a large magnetocristalline 
anisotropy, comparable in size to LIq type FePt films, were 



successfully fabricated ] i Consequently, we include also 
investigations of Lli ordered FePt and CoPt alloys in our 
study. 

By means of density functional theory (DFT) calculations 
we have determined surface energies and surface magnetism 
of various low-index surfaces, including the (100), (001), 
(110), (Oil), and the (111) facet in the LIq phase (cf. Fig. 13 
as well as the (111) surface of the Lli structure (cf. Fig.[T]). 
Regarding the surfaces of elemental systems, including bcc 
Fe surfaces^Oi^L^, Pt (1 1 1) and (001) surfaces^, and 4d 
transition metal surfaces^^i^ numerous studies can be found 
in literature , 38.39.40,4 1 p^j. i^jjjary transition metal alloys how- 
ever, only few investigations are available^^i^ 
To the best of our knowledge, this is the first systematic first 
principles comparison of the energies of various low-index 
surfaces of FePt, CoPt and MnPt with LIq and LI i order. 

n. METHOD 

For the evaluation of the surface energies, the so called slab 
approach^- was used. Here, the semi-infinite problem is rep- 
resented by a periodically repeated two-dimensional slab with 
two surfaces separating the periodic images by a sufficient 
amount of vacuum in the third direction. Some representative 
slabs are shown in Fig.|2l Here it can be seen that for binary 
alloys two different surfaces are encountered when slabs are 
stacked along the [001] and [110] direction in the LIq struc- 
ture and along [1 1 1] in the LI i phase. In those cases, one sur- 
face is entirely covered by Fe atoms while the other is covered 
with Pt atoms, respectively. Then the surface energy should 
be divided into two element-specific contributions as one sin- 
gle material component may be predominantly found at the 
surface. For these cases, surface-energy phase diagrams have 
been evaluated in order to account for the surface energies 
of the single material constituents. Regarding other metallic 
surfaces and, in specific, semiconductor surfaces, systematic 
investigations have been devoted to obtain structure, surface 
free energies, and segregation properties by, among others, 
the group of M. Scheffler4^iMiiM2i50^ 
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A. Computational details 
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FIG. 1: (color online) Left: The Llo-unit cell. Right: \A\- 
cell as used in the calculations. Dark (blue) spheres denote 
Fe/Co/Mn atoms, light (magenta) spheres Pt atoms. In the LIq or- 
der, monoatomic planes are stacked along the [001] direction, in the 
Lli structure, along the [111] direction. The tetragonal distorted 
LIq structure has two different lattice parameters a and c. In the LI i 
structure a slight distortion along the [111] direction may occur. The 
Lli crystal structure has only rhombohedral symmetry. 
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FIG. 2: (color online) Supercells (for clarity repeated in x and y - 
direction) used in the calculation of the surface energies in LIq FePt, 
CoPt and MnPt. As in Fig.[Tl dark (blue) atoms are Fe/Co/Mn, light 
(magenta) atoms are Pt. Top left: slab for the (001) and (110) sur- 
face; bottom: the (100), (Oil), and the (111) surfaces. Two different 
surfaces appear in [001] and the [110] directions. The two surfaces 
are covered by different elements, while in the other cases the slabs 
are limited by identical mixed surfaces. 



The self-consistent calculations are carried out with the 
Vienna ab initio simulation package (VASP) using a plane 
wave basis set and the projector augmented wave (PAW) 
framework] ^^-^^ The exchange-correlation potential is used 
in the functional form of Perdew, Burke and Ernzerhof 
(PBE)i-?^51 The PAW potentials include the following valence 
electrons: Fe: 3p^3d^4s', Co: 3dHs', Pt: Sd'^Ss' and Mn: 
3d^4si. All plane waves with energies below the cut-off en- 
ergy are included in the basis set. The cut-off energies were al- 
ways chosen 25% larger than the largest default cut-off of the 
element-specific potentials. We used (in eV): Fe: 366.5, Co: 
335.0, Mn: 337.3, FePt: 366.5, CoPt: 335.0 and MnPt: 337.3. 
The integration over the Brillouin zone is done by means of fi- 
nite temperature smearing (Methfessel-Paxton method for the 
surfaces) or tetrahedron method (for bulk systems). For the 
first case, the parameter a determines the width of the smear- 
ing (in eV): y-Fe: 0.15, a-Fe: 0.32, Co: 0.15, FePt: 0.2, 
CoPt: 0.28 and IVInPt: 0.2. We used the following k-point 
grids: For bulk calculations: A F centered (G) (13/13/13) 
grid. For (001) surface calculations: y-Fe: G (19/19/1), a-Fe: 
Monkhorst (M) generated (16/16/1) grid, fee Co: IVI (16/16/1), 
Pt: G (19/19/1), FePt: M (16/16/1), CoPt: IVI (14/14/1), and 
MnPt: M (14/14/1). The parameter o has been carefully cho- 
sen, so that the entropy term is lower than ImeV/atom. The 
electronic self-consistency iteration cycle is aborted when the 
energy difference between the old and the new energy is less 
than 10~' eV. The slabs used to model the surfaces consist of 
up to 32 atomic layers. Adjacent supercells are separated by 
a vacuum region of about 15 A to avoid interaction between 
neighboring supercells (Fig. |2l). The geometric relaxation is 
done by the conjugate gradient algorithm and at least the out- 
ermost 4 layers are optimized. Relaxation was stopped when 
the forces were less than 0. 1 eV/A. 



B. Surface energy calculations 

Surfaces can be created by dividing an infinite crystal into 
two parts. The energy needed to cut the bonds and bring the 
two resulting parts to infinity determines the surface energy. 
A straight forward procedure to calculate the surface energy 
is to examine the total energy E^°^(n) of a slab of the material 
of interest with n atomic layers and to subtract n times the 
bulk energy Ebuik of an atomic layer obtained from a separate 
calculation (e.g., see Ref. HUfor a detailed introduction): 



lim<! ^(£'°'(«)-«-£buik) 



(1) 



Here, A is the unit surface area. For sufficiently thick slabs, 
bulk properties are approached in the interior of the slab and 
y is expected to converge as a function of the slab thickness 
towards the exact surface energy. But proceeding as above, 
the surface energy rather diverges with slab thickness due 
to slight, unavoidable, numerical discrepancies which can be 
caused, e.g., by the choice of different basis or k-point grids in 



3 



slab and bulk calculationia^ In order to avoid this, the so called 
slab approach was proposedi^ Within the slab approach, the 
bulk energy, iibuik, is estimated from the same slab systems 
for the surface energy calculations instead of using a single 
separate bulk calculation within a small unit cell. Therefore, 
the divergence problem can be avoided as consistency in all 
technical parameters is maintained. 

The quantity /ibuik is extracted as follows: First all of the to- 
tal energy of various slabs with increasing thickness is cal- 
culated and plotted versus slab thickness. For large enough 
thicknesses the slope of a fitted straight line yields the bulk 
energy isbuik- In the present calculations, the surface energy 
converges properly if slabs with 6 atomic layers or less are 
discarded. 

When an Llo FePt crystal is cleaved on any of the (100), 
(Oil), and (111) planes, the two exposed surfaces are of mixed 
atomic composition, i.e. they consist of the same amount of 
Fe and Pt atoms. On the other hand, in the cases of both (001) 
and (110) cleavages, one surface consists entirely of Pt atoms 
while the other surface consists entirely of Fe atoms. Thus, 
proceeding as described above, we obtain an averaged value 
over the surface energies of both orientations but no informa- 
tion about the element specific contributions ype and ypt. How- 
ever, a variation range for the surface energies can be given by 
means of surface-energy phase diagrams . Two equiv- 
alent surfaces on top and on the bottom of the slab require to 
consider off-stoichiometric systems. If we compare energies 
of non-stoichiometric systems, the chemical potentials /j, of 
the single material constituents become involved: 

1{N,} = ^(e'"\N,]-Y^N,-h^. (2) 

Here, A^i is the number of atoms of the material component / 
and /Ji its chemical potential. Eq. ^ is considered here only 
at temperature T=0. At finite temperatures, the total energy 
has to be replaced by the Helmholtz surface free energy. A 
detailed thermodynamic derivation can be found in the liter- 
ature, e.g., see Refs. I45lf57l for more details. For the case of 
FePt, Eq. (|2|i reads 

Y(A^Fe,A^Pt) - ^ (£""(A^Fe,A^Pt) -A'FeA'Fe - A^Ptm) • O) 

The surface atoms are in equilibrium with the surrounding 
bulk reservoirs, which consist of the pure Fe or Pt metal and 
the underlying bulk alloy. Thus, the chemical potentials /jpe 
and /jpt are not independent, but related to the bulk alloy chem- 
ical potential A'FePt(buik) ~ 2 • iibuik, the bulk chemical poten- 
tials of the elemental constituents /Jp""' and /Jp""' ™d the heat 
of the alloy formation AZ/pePt: 

A'FePt(bulk) = A'Fe + m = A'Pe"' + A'p""' - A/ZpePt • (4) 

Unhke in their bulk equilibrium phases, the chemical po- 
tentials of the single material constituents within the alloy, 
A'Fe + A'Pt> are not known. However, one can eliminate one 
of them, e.g. /jpe: 

ypt = ^ {E'°\NF„Npt) - A^Fempt - AA^m) • (5) 
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TABLE I: Surfaee energies, a, in eV/atom and surfaee layer spin 
moment, M, in ^Ufi/atom of the faeets (111), (001), and (110) for a- 
and y-iron, fee cobalt and platinum. Unrelaxed (unr) and relaxed (r) 
geometries are compared. In the last column the bulk spin moment 
per atom is given. 

Here, we assume a slightly Pt-rich environment and the sur- 
face stoichiometry is given by ziN=Npt— Npe . The stability of 
the bulk alloy against decomposition requests that the chemi- 
cal potential /jpt can take only values in the range: 

A'Pt(bulk) ^ |A//FePt| < jfPt < A'Pt(bulk) (6) 

Now we can express Eq. (|5]l as a function of the difference in 
Pt chemical potential, jUpt ^A'pt(buik)- 

ypt = ^ (A^Fe , A^pt ) - A%e • feptibuik) - AA^A'pt(buik) 

-AA^[//Pt-A'pt(buik)])- (7) 

With the help of this equation, the so-called surface-energy 
phase diagrams can be determined. This approach has been 
applied successfully to estimate the stability of various semi- 
conductor surface reconstructions 4^i^2i^2i^ 



III. RESULTS 

A. Elemental Systems 

For the elementary systems, the properties of the bulk 
phases are well known and have been reported previously 
for most caseS] "^'i"^^'"^^i"^V^i"^V^i"^^ Our results match well with 
these investigations especially with those, where similar meth- 
ods and technical parameters were used i^^i^^'"^^'"^^ For met- 
als with fee lattice structure, we could confirm the follow- 
ing trend: With decreasing coordination number of the sur- 
face atoms, the surface energy, a, and the spin moment, M, of 
the outermost surface layer increases. This correlation is well 
known in literatur o^^i^^ and leads for fee metals to: 

o(lll) < o(OOl) < a(llO), 
M(lll) < M(OOl) < M(llO). 

These results match intuition since the surface energy (the en- 
ergy needed for cutting some "bonds") must grow with de- 
creasing surface coordination number The coordination num- 
bers (z) for the fee structure are: z~l for the most open sur- 
face (110), z = 8 for the (001) facet and z = 9 for the most 
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TABLE II: Summary of relaxed surface energies in eV/atom (J/m^) for fee Co, Pt, fee Fe, and bee Fe as calculated by the authors in compar- 
ison to data taken from the literature. 
° Authors: GGA (PBE), VASP 
* GGA (PW91), VASP 

Full charge density (FCD) method in the GGA, based on linear muffin tin orbitals in tight binding (TB-LMTO) and atomic-sphere approxi- 
mation (ASA). Unrelaxed surface energies, using a lattice constant a = 3.001A. 

LDA, FP-LMTO, and Ceperly-Alder parametrization for xc-potential, slab approach. 

VASP, TB-LMTO spinpolarized GGA for xc-potential. 
■f TB-LMTO, ASA, Ceperly-Alder for xc-potential. 

^ TB-LMTO, ASA, Ceperly-Alder for xc-potential. Surface relaxation neglected. Using experimental lattice constant. 

Tight binding method (TB). Unrelaxed surfaces geometries. 
' All-electron full-potential linearized augmented plane-wave FPLAPW (WIEN97), GGA(PBE). 

FPLAPW (WIEN97), GGA(PBE). Using the experimental lattice constant. 
^ FPLAPW, LDA, Barth and Hedin formula for exchange-correlation potential. 
' Empirical, modified embedded-atom method: MEAM. 

Extrapolation from experimental solid-vapour surface energies at higher temperatures to T = K (approximation for an "averaged (av)" 
polycrystalline surface). The solid-vapour surface energy is derived from liquid surface-tension measurements. 

" Estimation of surface energy by subtracting from the measured surface-tension of the liquid an entropy term propotional to the melting 
temperature. 



densely packed (111) surface orientation. For the very open 
bcc geometry we have: z = 4 for the (001) and the (111) sur- 
face and z = 6 for the (110) surface. Table |T] presents our 
calculation of surface energies, a, in eV/atom and magnetic 
moment on surface atom, M, in /j^/atom of all considered el- 
emental systems. Results for unrelaxed (unr) as well as re- 
laxed (r) structures are presented. In the last column, the bulk 
magnetic moment is given. For platinum, a large variation in 
surface energy (by a factor of two) between the most open 
(110) surface and the densely packed (111) surface is found. 
In the surface layers, the spin moment is enhanced between 
2 and 8% compared to the bulk value. This is in accordance 
with the observation of increased magnetic moments in low 
dimensional systems, e.g., as in small Fe clusters i ^^'^*^ 
In Table [III our results for the surface energies are compared 
to data available in the literature. For Pt, we find very good 
consistency of our surface energies with the ab initio calcu- 



lations of Refill [m Concerning the work of da Silva 
et al.^ the deviations remain in the range of 4-8%. They use 
the same functional form for the exchange correlation poten- 
tial (PBE) and their slabs are relaxed as well. Noteworthy 
deviations occur in comparison with Ref.'s?, where DFT cal- 
culations within the local density approximation (LDA) are 
performed. The authors use seven layer of vacuum in between 
adjacent supercells and do not relax their slabs. Their surface 
energy for the Pt (001) facet Opt(OOl) = 1.245 eV/atom is by 
27% larger than our value of apt(OOl) = 0.908 eV/atom, and 
presumably related to the different choice of the exchange- 
correlation potential. For 3c/-transition metals the LDA is 
known to show strong over-binding, i.e. cohesive energies 
turn out to be too large and lattice constants too small com- 
pared to experiment. As the surface energy is correlated to the 
cohesive energy, the values obtained within the LDA can be 
expected to be too large as well. This is different for the gra- 
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dient corrected exchange correlation potential where a slight 
overestimation of the lattice constant is common. 
The values reported by Skriver et alM using the tight- 
binding linear-muffin-tin orbital approach (TB-LMTO) with 
the atomic sphere approximation (ASA) are considerably 
larger than our PBE-values for all systems under investiga- 
tion, e.g., larger by 0.4 eV/atom for Co (111). Similar con- 
siderations hold true for the data of Ref.|65| which again show 
strong deviations to larger values. For completeness, we list 
in Table nil also semi-empirical and empirical methods, as e.g. 
the tight-binding and the modified embedded atom method 
(MEAM). These approaches do not yield the same accuracy 
of DFT methods but are frequently used for large scale sim- 
ulations. Here, free parameters are fitted to reproduce certain 
surface properties and thus can no longer be considered as 
high-level ab initio investigations. 

Also for a-Fe very good agreement with other first-principle 
DFT calculations (Ref.l6lll62[ Isol) is achieved. Slight devia- 
tions occur compared to the results of Ref.[63, 31, which 
may again be in part related to the different approximations for 
the exchange correlation potential as discussed above, differ- 
ent lattice constants or missing relaxation of the surface layer 
If experimental data are available at all, they mostly are ob- 
tained by liquid metal surface-tension measurement at higher 
temperatures and are then extrapolated to T = K. Thus, they 
are "averaged" values which can not be attributed to a special 
surface orientation. 

The large discrepancies found in literature show that the cal- 
culation of surface energies is a rather delicate task. They 
should be interpreted in terms of a comparison to other val- 
ues. This underlines the necessity for a systematic and (tech- 
nically) consistent comparative investigation of single element 
and binary transition metal surfaces as presented in this work. 



B. Binary Alloys of Fe, Co, Mn with Pt 

The structural and energetic properties of bulk Pt-based al- 
loys with Llo order have also been subject to numerous the- 
oretical and experimental surveysi ^i^i^''''^i^''^^'^^i^'^ However, 
as we also deal with the less well studied LI i structure, we 
provide a detailed comparison of both phases in the follow- 
ing. 



1. LIq and LI i bulk phases of FePt and CoPt 

In the left panel of Fig. [3] the energies of different magnetic 
structures of the ordered Llo and LI i phases are compared as 
a function of atomic volume. The ferromagnetic (FM) phase 
with Llo order is the ground state with a lattice constant of 
3.835 A. The layer- wise antiferromagnetic order (AF) is only 
13.6 meV/atom higher in energy. This shows the competing 
behavior between ferromagnetism and antiferromagnetism, 
which has been predicted previously from ab initio calcula- 
tions for Llo FePt.iiSiMiS^ The Lli structure is characterized 
by an equilibrium lattice parameter of a = 3. 844 A and is 
found to be 122 meV/atom higher in energy. This is in 



agreement with the experimental observation that bulk Lli 
FePt is not stable. However, the energy difference between 
the phases decreases as the valence electron concentration 
increasesi^ As expected, the total magnetic moment, Mtot, 
in the FM Llo structure and FM Lli structure is dominated 
by the Fe spin moment, Mpe, and steadily increases with 
increasing volume. The induced Pt moment, Mpt, follows 
the trend of the Fe moments. The optimum c/a-ratio is 
determined keeping the volume at the energetic minimum 
of the cubic structures fixed (cf. Fig. [3] right). For the 
ferromagnetic phase, c/a = 0.974 minimizes the total energy 
while the antiferromagnetic phase becomes stable at a slightly 
lower c/fl-ratio. For the Lli structure the b/a-rdl\o is varied 
and shows that b/a = 1.015 minimizes the total energy 
(cf.Fig.O. 

The total magnetic moment in the Llo structure decreases 
with increasing c/a-ratio, while in the Lli phase shows only 
little variation. Analogous bulk calculations have been also 
carried out for CoPt and MnPt. MnPt possesses an antifer- 
romagnetic groundstate with a = 3.887 A and c/a = 0.937. 
For bulk CoPt, the energy versus volume curve and the c/a 
-variation (and b/a -variation for the Lli structure) is shown 
in Fig.m We find the FM Llo structure to be the most stable 
one with an equilibrium lattice constant of a = 3.793 A and 
c/a ~ 0.976. In the FM Lli structure a lattice parameter of 
a = 3.801 A and a b/a-ratio of 1.017 is obtained. The FM 
LI 1 structure is 68 meV /atom higher in energy. The magnetic 
spin moments show qualitatively the same behaviour as for 
FePt. But interestingly, the induced Pt moment is as high as 
in the case of FePt, even though the spin moment of the Co 
atom is clearly lower than the spin moment of the Fe atom. 
Thus, the hybridization between the Co and the Pt J-electrons 
seems to be stronger than in the FePt alloy. Furthermore, for 
CoPt, the hybridization is stronger in the Llo phase than in 
the LI 1 phase. 



2. Surface properties of FePt, CoPt, and MnPt 

Most preceding surface energy studies for binary alloys 
have been carried out for ideal cleaved surfaces, neglecting 
the effects of possible relaxations. To close this gap and to 
give an account on the importance of relaxations in FePt and 
CoPt surfaces, the surface energies of Llo and Lli FePt and 
CoPt with subsequent relaxation of the atomic positions are 
compared in Fig.|5]for all investigated surfaces. For the slabs 
with two different surfaces the values are averaged over both 
possible terminations, i.e., for the Llo phase, the (001) and the 
(110) facet and in Lli the (111) facet. 

Lowest surface energies are found for the highly coordi- 
nated (111) facet: In the case of Lli FePt and CoPt we 
find: o(lll) = 0.675 eV/atom, in the case of Llo FePt: 
a(lll) = 0.701 eV/atom, and for Llo CoPt our calculations 
yield: a =^ 0.654eV/atom (cf. Table HIB. These results may 
to some extent explain the trend that in gas phase experiments 
frequently FePt icosahedral nanoparticles with platinum cov- 
ered (111) surfaces are generated.'Wi The more open (100) 
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FIG. 3: (color online) Left: Energy versus volume curves and element-specific magnetic moment of LIq FePt for different magnetic structures 
and for the LI i phase. The ferromagnetic phase (FM, black circles) leads to the equilibrium structure. The antiferromagnetic phase (AF, orange 
diamonds) appears only 13.6 meV/atom higher in energy and shows the competition between ferro- and antiferromagnetism. The ordered LI i 
structure (Llj, green squares), the MnPt-type antiferromagnetic structure (AF2, blue downward triangles), and the non-magnetic phase (NM, 
red crosses) are not stable for bulk FePt. Right: Energy (distortion) with fixed volume for LIq and LI i FePt and different magnetic states. The 
ferromagnetic phase minimizes the energy at a c/a-ratio of 0.974 while the antiferromagnetic structure becomes stable at a slightly smaller 
c/a-ratio. For the Llj structure the i/a-ratio is varied and the energy minimum is found al b/a = 1.015. In the upper panels the magnetic 
moments Mpe and Mpt (in /jb per atom) are given as well as their sum Mtot (in /jb per cell). 



facets, which occur in the Llo cuboctahedron (cf. Fig. \TT[ 
which is desired for magnetic data storage media, lie with 
a (GDI) ~ leV/atom higher in energy. The highest surface 
energies are found for the most open facets (Oil) and (110). 
Qualitatively similar results have been found for Llo PdZn 
and PtZn.= Again, the surface energy decreases with increas- 
ing coordination numbers analogous to the trend for the fee 
metals. The modification of the surface energy by relaxation 
is maximum for the most open surfaces (cf. Fig. |5]i. Here, 
a reduction of the surface energy of about 7% occurs while 
for the densely packed (111) facet the reduction amounts only 
to approximately 2%. The layer resolved atomic relaxation 
processes is shown in Fig. |6] The outermost layers of the 
slab move slightly inwards. Here, the relative displacement 
amounts to about 1 .25% for the Fe surface layer, about 2% for 
the Co surface layer, and about 2.3% for the Pt surface layers. 
In the interior of the slab, the relaxation shows an oscillat- 
ing behaviour which disappears for FePt beneath the fourth 



subsurface layer This agrees with previous findings for bcc 
Fe surfacesSS and FePt nanoparticlesS.. The relaxation in this 
(GDI) slab system is significantly lower than for the (111) facet 
of a Llo FePt cuboctahedron (ca. This agrees well 

with the experimental findings in the case of cuboctahedra, 
where no noteworthy relaxation of the (GGl) and (IIG) sur- 
faces is found. 

The relaxation behaviour of elementary metal surfaces 
has been subject to various studies in the past six 
decades . ^^i^^i^^i^*^!^''^^!^^ For transition metals, it has been at- 
tributed to the competing influence of the partial pressures 
arising from the localized cZ-bonds on the one side and the 
ip-electrons on the other, which are partially relieved at the 
surfaced However, the picture for the complete transition 
metal series is not uniform. While for most transition metal 
systems with a nearly half-filled (i-band strong inward relax- 
ation is observed, the effect diminishes towards the end of the 
series and eventually reverses sign for the 5d noble metals Pt 
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TABLE IIL Surface energies, a, of various low-indexed facets in eV/atom and magnetic moment on surface atom, M, in /iB^atom for FePt, 
CoPt, and MnPt. The c/a-ratios of the considered face centered tetragonal structures (fct) are: FePt: 0.974, CoPt: 0.976, MnPt: 0.937. Values 
are given for unrelaxed (unr) as well as relaxed (r) surfaces. For the (001) and the (110) surfaces of the LIq phase and the (111) facet of the 
Lli phase, only averaged (av) values over both terminations can be given. For the antiferromagnetic MnPt alloy the absolute value of the Mn 
atom, |MMn| is shown. For LIq FePt an excellent agreement with the results of Ref. IH is obtained (deviations of less than 3%). 
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FIG. 5: (color online) FePt (blue squares) and CoPt (green diamonds) 
surface energies for the (ill), (100), (001), (Oil), (110) facets in 
Llo, and the (111) facet in LI] structure. The corresponding coor- 
dination numbers, z, are given as well. Unrelaxed (open symbols) as 
well as relaxed values (filled symbols) are shown. For the special ori- 
entations with two different surfaces, averaged (av) values over both 
surface energies are given. This applies in LIq for (001) and (110), 
and in Lli for (111). The highly coordinated (111) surfaces in Lli 
phase and LIq structure are energetically clearly favored. The lines 
are only guides to the eye. 



and Au 



36,92,93 



The effective inwards relaxation, which we ob- 



serve for the binary FePt and CoPt surfaces is thus certainly 
influenced by the hybridization of the 3d and 5d electrons 
within the surface and subsurface layer. 

Now we turn to the question how the single material con- 
stituents contribute to the averaged values. The exact deter- 
mination of, for example, Opt(OOl) in FePt (cf. Fig.|2| is not 
possible within the slab approach due to the missing knowl- 
edge of the chemical potential of the material components in 
the alloy. But the range of variation for the surface energy can 
be given by means of surface-energy phase diagrams. Here the 
two limiting cases for apt(OOl) in FePt are calculated by arti- 
ficially varying the stoichiometry and with this the difference 
in Pt chemical potential A/jpt = /jpt — A'pt(buik) explained in 
section HTBI cf. Eq. (|5]l. 

Our results for the surfaces energies in FePt are collected 
in the surface-energy phase diagram shown in Fig. [T] The 
data are in perfect agreement with Hong et alM. (see Table 
Uni l which were also obtained using the simulation package 
VASP. In addition we consider facets with two different ter- 
minations and surfaces compatible with Lli order. Vertical 
lines in Fig. |7] and Fig. [8] mark the limiting cases for the dif- 
ference in Pt chemical potential A/jpt given by the formation 
enthalpies A/Zpept in the LI i and Llo phase and the zero value, 
as for negative values (A/j < 0) no alloying would occur. At 
the right side of the diagram, a Pt rich environment is assumed 
and thus /jpt = l^pt{buik) A/jpt = 0: The surface atoms are 
in equilibrium with the surrounding Pt-metal and the underly- 
ing FePt bulk reservoir. On the left border, at AT/pePt — —0.6 
eV, an Fe-rich environment in the Llo structure is assumed 
and therefore /Jpt — A'pt(buik) — ^FePt- The vertical line at 
A/jpt ~ —0.4 eV corresponds to Af/pept in the LI i phase. The 
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FIG. 6: (color online) Relaxed atomic positions of a 12-layer thick 
FePt and 14-layer thick CoPt (001) slab in terms of the relative dis- 
placement Ad/d(ooi)- The outermost layers move inwards by a few 
percent in all cases. After the fifth layer beneath the surface of FePt 
(sixth layer for CoPt), the oscillatory expansion and compression of 
the ideal bulk lattice parameter disappears. Open (blue) squares de- 
note Fe atoms, filled (green) diamonds Co atoms. Circles (magenta) 
belong to Pt atoms, in the FePt alloy as open symbols, in the CoPt 
alloy as filled symbols. 



two limiting cases for A/jpt (inserted into Eq. |7]l yield the cor- 
responding limiting values for the realistic value of Opt. When 
we artificially vary the stoichiometry by varying A/jpt on the 
horizontal axis, at the same time the difference in Fe chemical 
potential changes inversely (upper horizontal axis), as can be 
seen from equation (HJi. 

The surface energy of Pt covered, i.e. not mixed, (111) facets 
in the Lli structure can take extraordinarily low values (in 
eV/atom): 



0.19 < apt < 0.4. 



(8) 



These have direct consequences for FePt nanoparticle mor- 
phologies and confirm the results of previous ab initio cluster 
simulations;21 These predict radially onion-shell type ordered 
core-shell icosahedra to be energetically favored over single 
crystalline Llo cuboctahedra for small particle diameters. The 
low energy of Pt-terminated (111) facets overcompensates the 
high energy of the 20 twins, which possesses an individual 
Lli order and the additional contribution of the twin bound- 
aries. 

For comparison the particle surface energies of orientations 
with mixed atomic composition as the (100) and the (111) 
facet in the Llo phase are given in Fig. |7] too (horizontal 
lines). 

A further promising candidate for future ultra-high density 
magnetic storage devices is LIq CoPt due to its similar high 
magnetocristalline anisotropy energy in the bulk phase (Ku = 
5 • lO^erg/cm^). Therefore, analogous calculations were done 
for CoPt surfaces (cf. Fig. [HJ. Qualitatively the same trends 
for the surface energies were found. Pt-covered (lll)-facets 
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FIG. 7: (color online) Equilibrium surface-energy phase diagram of LIq and Llj FePt including all investigated surface orientations. On the 
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the case of iron on the upper horizontal axis, respectively). The maximum difference is given by the formation enthalpy, AHp 



-0.6305 



eV/atom in the LIq phase and AHp^p^ = —0.4 eV/atom in the Lli phase, and symbolized by the two vertical lines. Thick bright (magenta) 
lines denote Pt covered surfaces, corresponding thick dark (blue) ones Fe covered ones. In the LIq structure we use dotted lines for the most 
open (110) facet and dot-dashed lines for the (001) facet. The mixed (100) and (111) facets are shown as thin horizontal (black) lines. The 
(III) surfaces associated with the LI i bulk ordering are denoted by thick solid lines. 



in the Lli structure are here even more favorable (see Table 
Uni l. In general, we find that elemental, solely Pt-terminated 
surfaces are preferred over Fe covered and mixed surfaces of 
the same orientation. The consistently low surface energy of 
Pt covered facets may be regarded as one important driving 
force for the strong surface segregation tendency of Pt in these 
alloys i^^'^^ The lowest energy orientation for mixed surfaces 
is the highly coordinated (111) surface in accordance with the 
elemental systems discussed in Section lTlI Al 
For binary systems with perfect Llo-order, (111) surfaces 
which are covered by only one atomic species can not exist 
for geometric reasons. On the other hand, this surface modifi- 
cation can be realized for the LI i structure, which is however 
not stable for bulk FePt. Thus, a sufficiently low surface en- 
ergy may stabilize the LI i structure in small particles. 

In addition, we have investigated the distribution of the 
Fe (Co respectively) spin moments and induced Pt moments 
inside relaxed LIq and Lli FePt and CoPt slabs with var- 
ious surface orientations and different surface terminations 
(cf. Fig. |9]l. For the case of FePt, we consider an Fe termi- 
nation for the (001) slab in LIq order and the (111) slab in 
Lli order. (The (111) surface in the LIq order always con- 
sists of mixed atomic composition). For comparison we have 
chosen for CoPt a Pt-covered (001) and (110) slab in the Llo 



phase while the (111) slab in the Lli phase has one Co and 
one Pt surface. The (111) and (110) surfaces in LIq phase are 
qualitatively the same for both alloys and are thus shown only 
once. Again, the spin moment of the transition metal atom is 
enhanced by about 3 — 4% at the outermost layers. The in- 
duced Pt surface moments show a strong dependence of the 
number of neighboring transition metal atoms. In the outer- 
most layer, the Pt atoms lose a part of their magnetic partners. 
For the (111) FePt surface in the Llo phase, which consists 
of Fe and Pt atoms (mixed atomic composition), the missing 
magnetic Fe neighbors are decisive for the reduced Pt moment 
in the outermost layer. In contrast to the situation in the sub- 
surface layer: Here, the Pt atoms have full coordination and 
thus show an enhanced moment due to the large Fe moment 
in the surface layer. In a similar manner the Pt moment is 
slightly enhanced in the subsurface layer of the Fe terminated 
(001) slab in Llo order and the (111) slab in Lli order. For 
CoPt this simple rule does not seem to hold true as we find an 
enhanced Pt moment for the Pt surfaces in the (110) slab in 
Llo order as well. 

Recently performed ab initio cluster simulations revealed 
that the stability of single crystalline morphologies might 
be stabilized in these alloys by reducing the number of 3d 
electrons j^^i2i2^ In an extreme case, this may be achieved by 
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FIG. 8: (color online) Equilibrium surface-energy phase diagram of LIq and Lli CoPt for the (100), (001), (110), and the (111) surface 
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m 
3. 



3.05 

<D 3 

LL 

2.95 

0.45 
a 0.4 
FO.35 
0.3 



1 1 1 1 1 1 1 

o-o( 
y *-^if ( 

q \ v-^( 
" \ \ □-□ ( 

1 1 1 TT * 1 ' 


1 1 1 1 1 1 

001) /L1„ p 

iii)/Li, y ''^ 

111)/L1„ / p 

110) /Li„ / r 

n 

M * 1 T 1 1 


1 1 1 1 1 



1 1 1 1 1 1 

Co " 

1 T 1 1 1 1 


CD 




Pt 




T 


- 



1.95 ^cj 

1.9 

0.45 
0.4 

0.35 ^ 
0.3 



-8 -6 -4 -2 2 4 6 
Layer number 



-6 -4 -2 2 4 6 8 
Layer number 



FIG. 9: (color online) Distribution of spin moments in relaxed FePt and CoPt LIq and Lli slabs for various surface orientations. The (111) 
surface in the Lli phase is shown as stars (black), all other symbols refer to the LIq phase: Circles (blue) mark the (001) surface, triangles 
(red) the (111) surface, and squares (green) belong to the (110) facet in the LIq phase. For the (111) facet in the LIq order (mixed atomic 
composition), we used a seven layer slab which is plotted with doubled layer spacing (always one empty layer in between) for a more convenient 
comparison with the other surface orientations. Note the different scaling for the Co spin moments, Mco, on the right. 



11 



changing from FePt to MnPt, as Mn has one 3d electron less 
than Fe. On the other hand, one has to take care of the strong 
antiferromagnetic tendencies present in Mn alloys. Therefore, 
we considered the most relevant surfaces (001) and (111) in 
the Llo structure of AF MnPt. A collection of calculated 
surface energies and surface spin moments for FePt, CoPt 
and MnPt are listed in Table |III] The surface energies of 
the (001) facets only vary in between 15 meV/atom for the 
different binary alloys. For the (lll)-facet the variation 
amounts to 70 meV/atom. Also for MnPt the (111) surface is 
more favorable than the (OOl)-facet with an even increasing 
energy difference. Further investigation of non-stoichiometric 
ternary Fe-Mn-Pt alloys might thus be an interesting task. 



3. Stability range of LIq- versus Ll\-ordered clusters 

In order to get an idea of the influence of the calculated sur- 
face energies on the equilibrium shape of small nanoparticles, 
we apply a simple approach to approximate the stability range 
of different structural morphologies in the Llo and compet- 
ing LI 1 order. In the limit of large diameters the particles can 
be regarded as spherical. For a first rough estimate we make 
only use of the energy differences between the Llo and Lli 
order for volume and for surface atoms of the lowest energy 
surfaces. These are the platinum covered (111) facets in the 
LI 1 structure and the platinum covered (001) facets in the LIq 
phase. A more reahstic picture should also take into account 
twin boundary energies and internal stress. As the Lli phase 
is not stable for bulk FePt, completely Llo-ordered particles 
are expected for large diameters. But with decreasing volume, 
the surface to volume ratio increases and the extraordinarily 
low surface energy of the platinum covered (111) facet in the 
Lli phase gains increasing importance. Because of this, the 
Lli-order becomes more favorable than the Llo-ordering for 
particle sizes below a critical diameter. At this critical diame- 
ter, the gain in surface energy is equal to the energy loss due 
to LI 1 ordering. For FePt we find a critical diameter of 3.7 nm 
and 6 nm for CoPt for averaged values of o in agreement with 
total energy calculations of binary transition metal clusters ( 
AppendixlAlfor more details);^ 

However, the equilibrium crystal shape of an arbitrary particle 
is not necessarily spherical. Following the investigation of the 
structural stability of single crystalline and multiply twinned 
FePt nanoparticles, which has recently been performed by 
Miiller and AlbCfi-™? we apply a more detailed continuum 
model, in which the different surface energies of the various 
facets are taken into account in terms of a Wulff construction. 
Miiller and Albe considered symmetric particles and surfaces 
with mixed atomic composition only. We will also allow for 
asymmetric particle morphologies in the following, details are 
given in Appendix iBl Thus, an octahedron which is solely 
terminated by (111) facets may be considered as a favorable 
particle morphology (cf. Fig.llOb. Other candidates are mul- 
tiply twinned morphologies as icosahedra. Here, the calcula- 
tion of optimum shapes requires the calculation of twinning 
energies which is beyond the scope of this paper. Therefore, 




FIG. 10: (color online) Schematic view of a regular octahedron, ter- 
minated solely by eight (111) facets. 




FIG. 11: (color online) Final shape of the LIq ordered, asymmetric 
Wulff polyhedron for the case of FePt with two Pt covered square 
(001) surfaces (green) on top and bottom, four square (100) facets 
with mixed atomic composition, and six hexagonal (111) facets at 
the sides. The distances of the three different facets to the particle 
center (dooi , dioo, and dm) are determined following the Wulff con- 
struction making use of the calculated surface energies (see Table 
Ullt . The area of the Pt covered (001) surfaces is considerably en- 
larged compared to the mixed (100) and (010) surfaces. For CoPt the 
particle shape is qualitatively the same. 



we compare only two competing single crystalline structural 
motifs: On the one hand, the above mentioned Lli ordered 
octahedron with two elemental Pt covered (111) surfaces and 
six (1 1 1) facets of mixed atomic composition (cf. Fig.fTOland 
Fig.[T2li and on the other hand, the Llo ordered Wulff polyhe- 
dron with eight (111) facets of mixed atomic composition, two 
elemental Pt terminated (001) facets, and four (100) surfaces 
covered with Fe as well as Pt atoms (cf. Figure [TTli. 

The resulting energy differences, E^p — Eq^J^^, are shown in 




FIG. 12: (color online) Calculated shape of the Llj ordered, asym- 
metric FePt octahedron with two Pt covered, hexagonal (ill) Ft sur- 
faces (green) on top and on the bottom, and six hexagonal (ill) 
facets at the sides with mixed atomic composition (orange). The dis- 
tances of the two different ( 1 1 1 ) facets to the particle center (d^\ j and 
dill) are determined using the calculated surface energies (see Table 
llllt and applying the Wulff theorem. As in the Wulff polyhedron, the 
area of the Pt covered (111) surfaces is considerably enlarged due 
to their extraordinarily low surface energy. For the case of CoPt a 
similar shape is obtained. 



Figure [13] for FePt (blue) and CoPt (green). These are given 
as a function of c/m, describing the distance of a (1 1 1) facets 
from the particle center. In addition, we consider the possible 
variation of the surface energy of Pt terminated (001) facets in 
Llo and (111) facets in LI i order with the chemical potential. 
This it is important to mention, since these asymmetric binary 
structures are in general non-stoichiometric and the compo- 
sition differs between the two morphologies. Therefore, in a 
strict sense, we can only give an estimate of the stability range 
of the different morphologies in the two competing ordered 
phases. 

For FePt as well as CoPt the asymmetric L 1 1 ordered octa- 
hedron is the energetically preferred particle morphology for 
sufficiently small particle sizes. Using averaged surface en- 
ergies of Pt terminated (001) facets in Llo order and (111) 
facets in Lli order, the Lli octahedron is lower in energy 
for distances dm ~ 2.76 nm for FePt (blue solid line). This 
corresponds to particle diameters up to 6.28 nm (assuming 
a spherical particle with equal volume) and a total number 
of approximately 9200 atoms per cluster. The energetic ad- 
vantage of the LI 1 octahedron becomes maximum for di 1 1 ~ 
1 .8 nm which is equal to a particle diameter of approximately 
~ 4.12 nm (about 2600 atoms) and amounts to 93 eV/cluster 
(35.7 meV/atom) (see Figure [T3]l. Applying the same consid- 
erations to CoPt leads to a critical particle diameter of approx- 
imately 11.26 nm below which the Lli order is the thermo- 
dynamically stable phase. At a particle diameter of ~ 7.5 nm, 
the maximum energy difference (18.18 me V/atom) is reached. 
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FIG. 13: (color online) Difference in energy between the LIq or- 
dered, asymmetric Wulff polyhedron (shown in Fig.ll lb and the LI i 
ordered, asymmetric octahedron (shown in Fig. I12t as a function of 
the distance d\\\ for CoPt (green) and FePt (blue). The solid lines 
correspond to the assumption of an averaged surface energy for the Pt 
terminated (001) facets in LIq order and Pt terminated (111) facets in 
LI 1 order, the dashed lines belong to maximum and the dotted lines 
to minimum Pt surface energy values. In the positive regions between 
the dotted and the dashed lines, i.e. the dark (blue) hatched area 
for FePt and bright (green) hatched area for CoPt, the Lli octahe- 
dron is the energetically favorable morphology, for negative energy 
differences, the Lli Wulff polyhedron is more stable. In the cen- 
tral, double hatched (turquoise) region the FePt and the CoPt regions 
overlap. Taking averaged Pt surface energy values (solid lines), we 
find that for distances up to di 1 1 ~ 2.76 nm for FePt and dm ~ 5 nm 
for CoPt (this corresponds to effective particle diameters of 6.28 nm 
and 1 1 .26 nm assuming spherical particles with the respective atomic 
volume), the Lli ordered octahedron is lower in energy. Here, the 
energetic advantage of the Lli octahedron becomes maximum for 
dm ~ 1 .80 nm for FePt and dm ~ 3.35 nm for CoPt. At this sizes, 
the LI 1 octahedron is about 93 eV/cluster (35.7 me V/atom) for FePt 
and about 300 eV/cluster (18.2 me V/atom) for CoPt lower in energy 
than the LIq Wulff polyhedron. 



Here, the CoPt LI i octahedron contains about 16500 atoms. If 
one takes the minimum possible surface energies for Pt termi- 
nated (001) facets in Llo order and (111) facets in LI i order, 
the Lli octahedron is stable up to dm ~ 4.9 nm for FePt, 
which corresponds to a critical diameter of 11.18 nm, and 
dm ~ 7.3 nm for CoPt (critical diameter of 16.47 nm). Tak- 
ing the opposite case, i.e. maximal Pt surface energies, leads 
to the dotted lines. Here, the Llo ordered Wulff polyhedron is 
the thermodynamically stable morphology for particle diame- 
ters larger than dm ~ 1.5 nm (critical diameter 3.43 nm) for 
FePt and dm ~ 3 .5 nm (critical diameter 7 .9 nm) for CoPt. In 
summary, the LI i ordered, Pt terminated, asymmetric octahe- 
dron can be expected for particle diameters between 3.43 nm 
and 11.18 nm for FePt and between 7.9 nm and 1 6 .47 nm for 
CoPt. This large expectation ranges of ~ 7.5 nm for FePt and 
~ 8.5 nm for CoPt must be seen as a kind of error bar, due to 
the missing knowledge of the chemical potential of the single 
material components. 

Nevertheless, in comparison to the simplified considerations 
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above these results yield on average even larger diameters 
and thus support the prediction that LI i ordered nanoparticles 
with Pt covered (111) facets are a competitive particle mor- 
phology for small cluster sizes for both, FePt and CoPt binary 
alloys. Multiply twinned morphologies as icosahedra, which 
have been investigated in Ref.|94l might even further optimize 
the area of favorable Pt-covered (111) surfaces, while provid- 
ing a more spherical shape. This however is achieved at the 
expense of internal interfaces, which need to be considered 
separately. 



IV. CONCLUSION 

We have calculated surface energies and surface magnetism 
of various low indexed surfaces for the elemental systems 
bcc and fee Fe, fee Co, fee Pt and for the binary (fct) alloys 
FePt, CoPt, and MnPt. For Fe, Co and Pt we have considered 
the (001), (110) and the (111) surfaces. In addition, also 
the (100), (Oil) orientation in the LIq phase and the (111) 
surface in the Lli structure have been examined for the 
binary alloys. The surface energies were determined using 
the slab approach. For the special surface orientations in 
the binary alloys with more than one possible coverage, 
surface-energy phase diagrams have been evaluated in order 
to account for the surface energy contributions of the single 
material components. 

For all systems under investigation, (111) facets show the 
lowest surface energy. Especially Pt covered (111) surfaces, 
as found in Lli ordered FePt, possess an extraordinarily 
low surface energy which is considerably lower than the 
respective surface energy of pure Pt. This gives rise to the 
preferred appearance of Pt terminated core-shell icosahedral 
nanoparticles in gas-phase experiments and agrees well with 
the results of theoretical cluster calculations which show 
that platinum terminated, radially Llo-ordered core-shell 
icosahedra are energetically favorable .^212^ 
The surface energies of CoPt qualitatively follow the same 
trend as found for FePt: y(Lli/lll) < y(111) < Y(100) < 
Y(001) < y(011) < Y(110). For LIq MnPt, the surface 
energies of the (001) and the (111) facet lie in the same 
range as those of FePt and CoPt while (111) surfaces are still 
previleged. We may speculate that the addition of Mn to FePt 
should not substantially modify the relation of the surface 
energies. The surface energy of purely Pt-covered surfaces is 
always lower than the energy of the elemental surfaces of the 
corresponding 3d-metal. 

These results allow us to estimate the stability range of the 
most favorable particle morphologies in Llo and LI i order In 
a first simple approach, relying only on the energy difference 
between the Lli and the LIq structure for bulk and for 
surface atoms, we can derive a critical diameter below which 
the Lli phase may be stabilized. For FePt, this diameter 
is about 3.7 nm, which is in good agreement with ab initio 
cluster simulations.^ For CoPt we find a critical diameter of 
6nm. Similar crossover sizes were obtained within a refined 
continuum model, which allows to assess contributions of 
different faces more precisely. Assuming single crystalline 



particles, candidates for stable structures can be determined 
following the Wulff construction: Favorable single crystalline 
morphologies are the non-spherical Wulff polyhedron with 
Pt covered (001) facets for the Llo order and the asymmetric 
octahedron with Pt covered (111) facets in case of the Lli 
order. Comparing the total energy of those two structural 
motifs and using averaged Pt chemical potential yields a 
critical diameter of ~ 6 nm diameter for FePt and ~ 1 1 nm 
diameter for CoPt below which the LI i ordered, asymmetric 
octahedron is the energetically preferred structure. 
This underlines the central result of this study, that the 
extraordinarily low surface energy of elemental Pt terminated 
(111) facets of Lli bulk crystals may stabilize FePt and 
CoPt nanoparticles with LI i crystalline order for sufficiently 
small particle diameters - although the corresponding bulk 
materials are unstable in the Lli structure. 
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APPENDIX A: SPHERICAL MODEL FOR THE CLUSTERS 

In the limit of large diameters the particles are assumed to 
be spherical and the difference in surface energy per atom is 



approximated by AEs = Op' "(111) — Op^ ' (00 1 ) . The volume 
energy difference per atom between the L 1 1 phase and the 



Llo phase is AEy ~ 



LIq 

bulk ' 



.LI, 
"bulk- 



The gain in surface energy 



is equal to the energy loss due to LI i ordering if the following 
condition is fuUfiUed: 



AE =N-AEw + S-AEs = 0. 



(Al) 



Here is the total number of atoms in the particle and S the 
number of surface atoms. Relation dAlb gives the percentage 
of surface atoms that leads to the stability of the Lli phase: 



AE^ 



S 
N 



(A2) 



With the bulk and surface energy differences taken from sec- 
tion|IIIB2J AEy = -0.13eV/atom and AEs = 0.39eV/atom 
for FePt, we estimate that the critical percentage of surface 
atoms per particle (8s) amounts to 32.7%. We can find the 
corresponding total number of atoms per particle using 
5s = 4/^.'^^ For 5s = 0.327 this leads to N~ 1840. The 
particle diameter can then be estimated with the help of the 
averaged atomic volume in Llo FePt, i2, using the formula 



APPENDIX B: ASYMMETRIC CONTINUUM MODEL 

As the surface energies of the various facets in LIq- and 
LI 1 -ordered FePt and CoPt differ considerably, the assump- 
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tion of a spherical particle is not necessarily valid. Rather 
asymmetric particles with an enlarged area ot the energetically 
favored facets are expected. This is taken care of in a more 
detailed continuum model. For single crystalline metal parti- 
cles, the thermodynamically stable shape is determined by a 
Wulff construction^^, where the energy minimizing shape is 
given by a constant ratio Yhu/t/hki with lihki the distance from 
particle center of a (hkl) facet with surface energy per unit 
area Yhki. The Wulff theorem applies to a macroscopic crys- 
tal. The ratios Tioo = Yioo/Yiii and Tooi = 1^oi/Yiii for the 
Llo structure, and f = Y^'i i /Yni for the LI i phase determine 
which particle morphology possesses the energy mini mizing 
shape. In the strong faceting limit (Yiio/Yih ^ ■\/3/2) and if 
the condition V3/2 < Fioo < is fulfilled, the Wulff shape 
is a truncated octahedron terminated by ( 111 ) and ( 1 00) facets 
only (cf. Fig.fTTTi. This applies to Llo FePt and CoPt particles 
where Fioo = 1.205 and Fioo = 1.265, respectively. In the 
limiting case of Fioo = — 1-732 the Wulff construction 
leads to a regular octahedron as depicted in Fig. [10] The Llo 
ordered Wulff polyhedron can be constructed with the help of 
the ratios F i oo = Yi oo /Yi 1 1 and Foo i = Yoq i /Yi 1 1 by trunc ating 
the vertices of a regular octahedron at distances dm, dioo and 
dooi from the centeri^^iiffi The volume V of the truncated oc- 
tahedron can easily be derived by subtracting from the total 
octahedron volume, Vocta, half of the volume of two small oc- 
tahedra truncated at the vertices in [001] direction, Vgcta' 
four small octahedra truncated in [100] directions, Vg^ta ■ 



V = Voc 



V 



dooi 
Octa ' 



(Bl) 



Aooi(diii) = 6d 



111 



Aioo(diii) =6dfii ( 1 - ^ 



(B6) 



(B7) 



Aiii(diii) = -\/3d^ii 



^^^"■[■-('-Tf)-K'-^) 

(B8) 

As the surface energy of Ft terminated (001) facets, t^qj, 
is a function of the difference in Pt chemical poten- 
tials, /jpt — A'pt(buik)' (Fig- 111 and Fig. [8]) the value of 
Fooi = "^^qj/yi 11 also varies between two limiting cases: 



7^0 



001 



mm 
001 

Ylll 



< 1 001 < i 001 = 



max 
001 



Ylll 



(B9) 



leading 
0.816 < 



to 0-686 < Fooi < 
Fooi < 0-950 for CoPt- 



1.074 for FePt and 
The averaged values 



are almost the same for FePt and CoPt: Fqqj = 0.88 and 
Fqq[ = 0.883, respectively. Thus, also the volume, V, and 
the areas Aooi and Am are not exactly determined but vary 
as a function of /jpt — A'pt(buik) - For the quantitative examples 
and representative particle shapes shown in section 1111 B 31 
averaged values for Fooi are used- 

The non-spherical, asymmetric LI i -ordered octahedron is 
defined by the ratios 



with Vbcta 



:4V3d^ii andy°™ 



di, 
Octa 



:4V3(d 



111 



ogously for V^ti) 



d|00\3 



(anal- 



The area of the (100) and (001) surfaces 
is simply given by the square of the edge length while the 
remaining (111) surface areas of the octahedron can be calcu- 
lated following the same idea when determining the volume: 



1 



Mil 



Ao, 



'cta (Ao^a - 



o AdlOO 
' ^^Octa 



)) 



with 



and 



Aocta = 12\/3d]ii 



jdooi 
*Octa 



12\/3(d 



111 



dooi n2 

' 



(B2) 



(B3) 



(B4) 



(correspondingly for Aq™^)- 
With this preliminary considerations and after substituting 
dioo = Fioodm and dooi = Fooidm, the volume V and the 
total area Ahki of (1 1 1), (100) and Pt covered (001) surfaces of 
the Llo Wulff polyhedron can be written as a function of the 
variable dm: 



V(dm)=4d?iiV3 



Fooi Y 



2 1 



i 100 \ 



(B5) 



7^1 



HI 



and 



Ik 
Tin 



0.444 for FePt 



0.446 for CoPt. 



(BIO) 



(BID 



As for Pt terminated (001) facets in Llo order, here, the Pt 
covered (111) facets are known only in between two limiting 
cases (Fig. ID and Fig. Ell yielding: 0.277 < f < 0.590 for 
FePt and 0.338 < F < 0.503 for CoPt. In analogy to the 
Llo ordered Wulff polyhedron, we express the volume V, the 
total area of Pt covered (111) facets A^'Ji, and (111) facets of 
mixed composition A™]' for the Lli ordered octahedron as a 
function of dm and F by 



V(dm): 



F 9 



AS'?(diii) = ^V3d?Hf, 



Ani(diii) = 4V3d2ij(3-r 



(B12) 



(B13) 



(B14) 



Under the assumption of an averaged value of the surface en- 
ergy, Y^lr = 0.733 for FePt and Y^'lt" = 0-722 for CoPt we 
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were able to precise the structural motifs. The resulting parti- 
cle shapes are shown in Fig. [TT] and [12] respectively. 
As expected, the Pt covered (001) facets are considerably en- 
larged by almost a factor of Aqoi / Aioo = 2.6 compared to the 
(100) facets with mixed atomic composition (Aqoi/Aioo = 3.3 
for CoPt). The ratio of the distance of the Pt terminated (001) 
facet and the mixed (100) facet from particle center (dooi 
and dioo, respectively) giving the aspect ratio of the particle, 
amounts to dooi/dioo = 0.73 (dooi/dioo = 0.70 for CoPt). 
Using the continuum model, also the total energy of a parti- 
cle can be expressed as a function of the distance of a (111) 
facet from the particle center dm. For large enough particles, 
it can be approximated by the sum of volume and surface en- 
ergy terms.i^ If denotes the total number of atoms in a 
particle, the particle volume is given by V = NQ., where Q. is 
the atomic volume. Thus, the total number of atoms varies 
with cluster size as 

N{dni)=V{dni)/^. (B15) 

If we neglect the twin boundary energy and the contributions 
of edge and corner atoms, we obtain: 

E{dui) ^N{dni)Ebu\k + Y,'^hki(.din)yhki- (B16) 

hkl 

Applying Eq. (IBI6I 1 yields for the Llo Wulff polyhedron 
(WP) 

£^'p« =NE^X + ^Anani +2Aooi7ooi +4AiooYioo> (B17) 



The dependence on dm is formally obmitted for simplicity. 
As mentioned above, also the total energy is a function of the 
difference in Pt chemical potentials, /jpt — A'pt(buik)' is de- 
termined only in a certain range. 

Analogously such considerations also apply for the LI 1 or- 
dered octahedron. Again, the low surface energy of Pt covered 
(111) facets will lead to an enlargement, compared to (111) 
facets with mixed atomic composition. We here also consider 
a non-spherical, asymmetric shape. Applying Eq. ( IB 16b to 
the LI 1 ordered, asymmetric octahedron gives 



EolL = NE^ui + 2Aruyl'n + 6A 1 n Yi 1 1 • (B 1 8) 



Under the assumption of an averaged value of the surface 
energy for FePt and CoPt we are able to predict structural 
motifs. Their shape is shown in Fig. [TJ] Indeed, the 
hexagonal Pt terminated (111) facets on top and bottom 
are enlarged by a factor of A^'j/A'™ = 3.44 for FePt 
(Afii/A™^ = 3.35 for CoPt) compared to the (111) facets 
of mixed atomic composition on the side of the particle. 
The distance of the Pt terminated (001) facet from particle 
center, dooi, is even more shortened compared to the mixed 
(100) facet, dioo, as found for the Llo WP. The ratio amounts 
to dfn/diii = 0.412 for FePt and dfjj/dm = 0.421 for CoPt. 
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